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Abstract. The oxyge n and nitrogen abundances in the Hll regions of the nine Virgo spirals of the sample from 


Skillman et al. 199t and in nine field spiral galaxies are re-determined with the recently suggested P - method. We 


confirm that there is an abundance segregation in the sample of Virgo spirals in the sense that the H l deficient 
Virgo spirals near the core of the cluster have higher oxygen abundances in comparison to the spirals at the 
periphery of the Virgo cluster. At the same time both the Virgo periphery and core spirals have counterparts 
among field spirals. Some field spirals have Hi to optical radius ratios, similar to that in Hi deficient Virgo core 
spirals. We conclude that if there is a difference in the abundance properties of the Virgo and field spirals, this 
difference appears to be small and masked by the observational errors. 
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1. Introduction 


There is evidence that 
properties of galaxies 


the environment affects 
in clusters (see review 


the 
of 

Balkowski 1992 ). The most obvious effect is probably 
t hat concerning the Hi c o ntent of galaxies in clusters 
(Haynes fc Giovanelli 1986, 


Huchtmeier & Richter 1989 


Cayatte et al. 1994 ). It is well established that the spi¬ 
ral galaxies of the Virgo cluster have a tendency to be H i 
deficient in comparison with normal field spirals (Solanes 
et al. 1996, 2001) and this deficiency is correlated with 
distance to the cluster center, the proportion of gas-poor 
spirals increasing continuously towards the cluster center. 
Thus, a cluster galaxy (especially a galaxy near the center 
of the cluster) evolves in a different surrounding gaseous 
medium in comparison to a field galaxy. 

The gas exchange between a galaxy and its ambient 
medium (loss of gas by the galaxy or gas infall onto the 
galaxy) changes the course of the chemical evolution of 
the galaxy. Then, it can be expected that the environ¬ 
ment affects the chemical properties of galaxies in clusters. 
Oxygen plays a key role in understanding the (chemical) 
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evolution of galaxies. The origin of oxygen seems to be re¬ 
liably established in contrast to other elements like carbon 
or nitrogen. The oxygen abundance can be considered as 
a tool to investigate the evolution of galaxies. For exam¬ 
ple, the value of oxygen abundance in a galaxy combined 
with the value of the gas mass fraction can tell us about 
the efficiency of mass exchange between a galaxy and its 
environment (Pilyugin & Ferrini 1998, 2000). 


A number of works have been devoted to searching for 
the possible effects of cluster environment on the chemi¬ 
cal properties of spiral galaxies (Shields et al. 1991; Henry 
et al. 1992, 1994, 1996; pkillman et al. 1996 ). The conclu¬ 
sions of these works are: i) the spirals at the periphery 
of the cluster are indistinguishable from the held galax¬ 
ies, ii) the Hi dehcient Virgo core galaxies have larger 
oxygen abundances at a predetermined galactocentric dis¬ 
tance r = 0.4i?25 (where R 25 is the isophotal radius) than 
the field galaxies of comparable luminosity or Hubble type. 


Accurate oxygen abundances can be derived from 
measurement of temperature-sensitive line ratios, such 
as [OHI]4959,5007/[OIII]4363. This method will be re¬ 
ferred to as the Tg - method. Unfortunately, in oxygen- 
rich Hii regions the temperature-sensitive lines such 
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as [OIII]4363 are too weak to be detected. For such 
Hii regions, empirical abundance indicators based on 
more readily observable lines were some years ago sug¬ 
gested ( Pagel et al. 1979 , Alloin et al. 1979 ). The empir¬ 
ical oxygen abundance indicator R 23 = ([OII]3727,3729 
-I- [OIII]4959,5007)/H^, suggested by Pagel et al. 1979, 
has found widespread acceptance and it has been used 
for the oxygen abundance determination in H ii regions 
where the temperature-sensitive lines are undetectable. 
This method will be referred to as the R 23 - method. 
Using the R 23 - method, the characteristic oxygen abun¬ 
dances (the oxygen abundance at a predetermined galac- 
tocentric distance) and radial oxygen abundance gradi¬ 
ents were obtained for a large samp le of field spiral galax¬ 
ies ( Vila-Costas & Edmunds 1992, Zaritsky et al. 1994 , 
van Zee et al. 1998 among others). 


abundance indicator is open to discussion (Zaritsky 1992, 
Kinkel & Rosa 1991, among others). It has been found 
(Pilyugin 2000) that the error in the oxygen abundance 
derived with the R 23 - method involves two parts: the 
first is a random error and the second is a systematic error 
depending on the excitation parameter. A new way of oxy¬ 
gen abundance determination in H ii regions (P - method) 
has been recently suggested (Pilyugin 2000, 2001a). By 
comparing oxygen abundances derived through the - 
method in high-metallicity H ii regions with those derived 
through the P ~ method, it has been found that the pre¬ 
cision of oxygen abundance determination with the P - 
method is comparable to that of the Tg - method (Pilyugin 
2001a,b). It has been also shown that the R 23 - method 
provides more or less realistic oxygen abundances in high- 
excitation H II regions but it produces overestimated oxy¬ 
gen abundances in low-excitation Hii regions. Taking into 
accou nt this fact t o gether with the fact known for a long 
time ( ^earle 1971 , Smith 1975 ) that galaxies can show 
strong radial excitation gradients, in the sense that only 
the low-excitation Hii regions populate the central parts 
of some galaxies, one can expect that the oxygen abun¬ 
dances of the inner Hii regions and the gradient slopes 
based on the ( 0 /H)fl 23 data can be appreciably overesti¬ 
mated. This speculation has been confirmed by compar¬ 
ison of the radial ( 0 /H)/j 23 abundance distribution with 
the radial (O/H)^^ abundance distribution within the disk 
of the well-observed spiral galaxy MlOl ( Pilyugin 200 Ibl ). 

The problem whether the cluster environment affects 
the chemical evolution of galaxies will be considered here 
based on the (0/H)p and (N/H)p abundances obtained 
in Section 2. In §3 we analyze these data with the mul¬ 
tiphase model applied to three Virgo galaxies, considered 
as typical examples. A dicussion is included in Section 4 
and Conclusions are in Section 5. 


2. Virgo and field spiral abundances 


2.1. Oxygen 

The (0/H)p oxygen abundances in Hii r egions are deter - 
mined with the expression suggested in (Pilyugin 2001a) 


12+log{0/H)p = 


R 23 + 54.2 -fi 59.45P -f 7.31P^ 


( 1 ) 


6.07 -b 6.71P -b 0.37P2 + 0.243i?23 

where R23 =i ?2 + Psj R 2 = I[OII]\372T+X3729 / I3, R3 = 
IIOIII]X4959+X5007/Ih/3, and P = i?3/i?23- 

For comparison, the (0/H)p23 oxygen abundances 
are also determined. Several workers have suggested 
calibrations of the R 23 in terms of the oxygen 
abundance (Edmunds fc Pagel 1984 , McCall et al. 1985 


However, the basic problem whether R 23 is an accurate 


Dopita & Evans 1986, Zaritsky et al. 1994, among oth¬ 
ers). Abundances for a large sample of spiral galaxies 
were derived by Zaritsky et al. (1994). The Zaritsky 
et al’s calibration is an average of the three calibra¬ 
tions by Edmunds & Pagel 1984, McCall et al. lOS"^ , and 
Dopita & Evans 1986. The oxygen abundances in Virgo 
spiral galaxies were determined by Skillman et al. (1996) 
using this same empirical R 23 - abundance calibration. 
Then this latter R 23 calibration has been adopted here for 
the determination of the ( 0 /H)p 23 oxygen abundances. 

The set of available spectra from the literature has 
been used for the determination of the (0/H)p and 
(0 /H)p 23 oxygen abundances. There exist nine Virgo spi¬ 
rals with line intensity measu rements for at least four H ii 
regions (^killman et al. 1996 ). We will compare the abun¬ 
dance properties of this Virgo sample to a sample of field 
spiral galaxies. This sample of field spiral galaxies also con¬ 
sists of nine galaxies, which were chosen in such way that 
1) they lie in the same range of luminosity as the Virgo 
sample spirals, 2) they cover the same range of morpho¬ 
logical type as the Virgo spirals, and 3) they have large 
enough number of H ii regions with measurements of oxy¬ 
gen and nitrogen lines. The adopted and computed pa¬ 
rameters of Virgo and field spiral galaxies are summarized 
in Table |^. The NGC number is listed in column 1. The 
adopted distance is reported in column 2 (source for the 
distance in column 3). The blue luminosity of the galaxy is 
reported in column 4. The isophotal radius R 25 in arcmin 
and the n umerical Hubble type of th e galaxy (T type) 
taken from de Vaucouleurs et al. 1991 (RC3) are given in 
columns 5 and 6 , correspondingly. The central (0/H)p 
oxygen abundance and the gradient expressed in terms of 
dex/R 25 are listed in columns 7 and 8 . The source(s) for 
the line intensity measurements is reported in column 9. 
The central (0/H)p23 oxygen abundance and the gradi¬ 
ent expressed in terms of dex/R 25 are listed in columns 
10 and 11 , respectively. 

Fig.|^ shows the derived (0/H)p and (0/H)p23 oxygen 
abundances and excitation parameter P for H ii regions 
in Virgo spirals as a function of galactocentric distance 
(normalized to the isophotal radius). The radial (0/H)p23 
gradients in the Virgo spiral galaxies derived here are 
close to those determined by Skillman et al. (1996), with 
the exception of the galaxy NGC 4651. Eor this galaxy, 
those authors have found a value of the gradient of -0.42 
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Table 1. The galaxy sample. 




P method 


R 23 method 

Galaxy 

d ref 

logLs R 25 T type O/Ho gradient 

literature sources 

O/Ho gradient 


Mpc 

(’) dex/R25 

for the spectra 

dex/R25 







Virgo cluster galaxies 




NGG 4254 

16.8 

SKSZ 

10.61 

2.81 

5 

8.94 

-0.65 

HPC,MRS,SSK 

9.44 

-0.65 

NGG 4303 

16.8 

SKSZ 

10.62 

3.23 

4 

8.84 

-0.72 

HPLC,SSK 

9.43 

-1.05 

NGG 4321 

16.8 

SKSZ 

10.66 

3.79 

4 

8.86 

-0.37 

MRS,SSK 

9.36 

-0.38 

NGG 4501 

16.8 

SKSZ 

10.70 

3.54 

3 

8.99 

-0.52 

SKSZ 

9.48 

-0.45 

NGG 4571 

16.8 

SKSZ 

9.94 

1.86 

6 

8.90 

- 0.20 

SKSZ,SSK 

9.34 

- 0.10 

NGG 4651 

16.8 

SKSZ 

10.24 

2.04 

5 

8.71 

-0.64 

SKSZ 

9.25 

-0.85 

NGG 4654 

16.8 

SKSZ 

10.34 

2.45 

6 

8.84 

-0.77 

SKSZ 

9.34 

-0.84 

NGG 4689 

16.8 

SKSZ 

10.10 

2.18 

4 

8.88 

-0.41 

SKSZ 

9.40 

-0.39 

NGG 4713 

16.8 

SKSZ 

9.88 

1.35 

7 

8.71 

-0.73 

SKSZ 

9.24 

- 1.02 






Field 

galaxies 





NGG 628 

9.7 

ZKH 

10.25 

5.36 

5 

8.71 

-0.48 

MRS, vZ, FGW 

9.31 

-0.77 

NGG 1232 

21.5 

vZ 

10.69 

3.71 

5 

8.73 

-0.58 

vZ 

9.36 

- 0.88 

NGG 2903 

8.9 

DK 

10.45 

6.29 

4 

8.94 

-0.73 

MRS, ZKH, vZ 

9.44 

-0.76 

NGG 3184 

8.7 

vZ 

9.93 

3.71 

6 

8.97 

-0.63 

MRS, ZKH, vZ 

9.53 

-0.80 

NGG 3351 

8.1 

ZKH 

9.90 

3.79 

3 

8.94 

-0.31 

MRS, OK, FGW 

9.43 

-0.33 

NGG 5194 

7.7 

ZKH 

10.49 

5.61 

4 

8.92 

-0.40 

MRS, D, BKG 

9.43 

-0.45 

NGG 5236 

4.5 

S 

10.31 

6.59 

5 

8.75 

-0.24 

DTJS, WS 

9.26 

-0.27 

NGG 5457 

7.5 

KG,vZ 

10.67 

14.42 

6 

8.76 

-0.77 

KG 

9.37 

-1.52 

NGG 6946 

5.9 

KSH 

10.60 

8.30 

6 

8.73 

-0.49 

MRS, FGW 

9.24 

-0.62 


List of r eferences'. 
BKG - 
FGW - 


KSH - 
SKSZ - 


Bresolii^_e£_aL_^99£ : D 
jFergurson et^LJh98 


DTJS - Dufour et al. 1980 


HPG 

iKarachentsev et al. 200C; MRS - 




Dfazet_ah_1991 ; DK - Drozdnvsk^^_J£arachentsev_2000 

HenrT'etal. 199 j; HP LC - ; KG - iKennic utt fc Garnett 1996 


_ killman et al. 

Zaritsky et al. 1994| ; 


1996 


SSK - shields et al. 1991 


McCall et al. 1985 ; OK - Oey fc Kennicutt 19931; S - schnijdt_etjik 


vZ - van Zee et al. 1998 


WS - [Webster fc Smith 1985; 


1994 


ZKH - 


(dex/R 25 ), while a value -0.85 (dex/R 25 ) is obtained here. 
The Hii region +131,+021 with galactocentric distance 
I.O 8 R 25 is responsible for the small value of the gradient 
derived by those authors for that galaxy. Close examina¬ 
tion of the radial distribution of (0/H)p oxygen abun¬ 
dance within the disk of NGC 4651 (Fig.|^ shows that 
the oxygen abundance at galactocentric distances larger 
than ^ O. 8 R 25 is expected to be below 12+log(0/H) = 
8.2. Therefore it is unlikely that the Hii regions within 
the disk of NGC 4651 at galactocentric distances larger 
than ^ O. 8 R 25 belong to the upper branch of the 0/H - 
R 23 diagram and, as consequence, the R 23 calibrations of 
the upper branch of the 0/H - R 23 diagram should not 
be used for the oxygen abundance determination in the 
Hii region +131,+021. Therefore this region has been ex¬ 
cluded from consideration here, resulting in a larger value 
of the gradient as compared to the previous value of the 
gradient in Skillman et al. (1996). 

As we can see in Fig.|^ most of the Hii regions ob¬ 
served in the Virgo spiral galaxies are of low-excitation 
and they do not show (or show marginal) radial excitation 
gradients. Therefore the (0/H )p gradients are close to the 


(0/H)p 23 gradients, but the abundance (0/H)p at any 
given galactocentric distance is significantly lower than 
the value estimated from the (0/H)p23 calibration. This 
reflects the fact that the R 23 - method may provide more 
or less realistic oxygen abundances for high-excitation H ii 
regions but produce overestimated oxygen abundances in 
low-excitation H ii regions. As an example, we see that a 
more or less distinct radial excitation gradient can be seen 
in the Virgo spiral NGG 4303. In consequence the (0/H)p 
gradient in the NGC 4303 is steeper than the (0/H)p23 
gradient, Fig.0. 

Fig.|| shows the derived (0/H)p and (0/H)p23 oxygen 
abundances and excitation parameter P for H ii regions in 
field spiral galaxies. Some of them (NGC 628, NGC 1232, 
NGC5457) show an appreciable radial excitation gradi¬ 
ents and, as a consequence, the (0/H)p gradients in these 
galaxies are steeper than the (0/H)p23 gradients, FigJ|. 

Fig.^ shows oxygen abundances at r=0.4 R 25 as a func¬ 
tion of blue luminosity (a) and as a function of the mor¬ 
phological type T (b) for the Virgo and field spirals. The 
oxygen abundances at r=0.4 R 25 have been used instead 
of the central oxygen abundances because the central oxy- 















































4 


L.S. Pilyugin et al.: Oxygen and nitrogen abundances in Virgo and field spirals 




Fig. 3. Oxygen abundance at r=0.4 R 25 as a function 
of blue luminosity (a) and as a function of T type (b) 
for the Virgo and field spirals. The pluses correspond 
to the field galaxies (NGC 628, NGC 1232, NGC 2903, 
NGG 3184, NGG 3351, NGG 5194, NGG 5236, NGG 5457, 
NGG 6946), the squares represent the Hi deficient (Virgo 
core) spirals (NGG 4501, NGG 4571, NGG 4689), the tri¬ 
angles represent the intermediate galaxies (NGG 4254, 
NGG 4321, NGG 4654), and the open circles represent 
the Virgo galaxies with normal Hi disks (NGC 4303, 
NGG 4651, NGG 4713). 


gen abundances in these galaxies can be affected by large 
uncertainties: the number of Hii regions in the Virgo core 
galaxies is low and they cover a small interval of galacto- 
centric distances, see FigJl, implying that the extrapola¬ 
tion method used to estimate these central data is uncer¬ 
tain. Since the value of r=0.4R25 is within the interval of 
galactocentric distances covered by the observed Hii re¬ 
gions in each galaxy, we believe the oxygen abundance at 


Fig. 4. 12-|-logO/H at r=0.4R25 versus Rh/R 25 diagram 
for Virgo (circles) and field (pluses) spiral galaxies. 

r=0.4R25 is more reliable than the central oxygen abun¬ 
dance. 

Skillman et al. (1996) divided the nine Virgo spi¬ 
rals in their sample into three groups of three galax¬ 
ies each: those with strong Hi deficiencies (NGG4501, 
NGC4571, and NGG4689), intermediate cases (NGC4254, 
NGC4321, and NGG4654), and those with no Hi deficien¬ 
cies (NGC4303, NGC4651, and NGG4713). The three sub¬ 
groups into which the sample has been divided are coded 
with different symbols in Fig.|[ The pluses correspond to 
the field galaxies, the squares represent the Hi deficient 
(Virgo core) spirals, the triangles correspond to the inter¬ 
mediate galaxies, and the open circles represent the Virgo 
galaxies with normal H i disks. 

The examination of Figs.^,b shows that there is an 
abundance segregation in the sample of Virgo spirals in the 
sense that the H i deficient Virgo core galaxies are all near 
the top of the abundance distributions of galaxies with 
similar properties, such as it was obtained by Skillman 
et al. (1996). At the same time, some field spiral galaxies 
have oxygen abundances similar to those ones from the 
most oxygen-rich (H i deficient) Virgo spirals. 

Those authors also found that the characteristic oxy¬ 
gen abundance (oxygen abundance at r=0.4 R 25 ) for the 
Virgo cluster galaxies correlates with the Hi to opti¬ 
cal radius ratio, Rir/R 25 , where Kh is the Hi radius at 
the isophotal level of 1 Mq pc“^ after Warmels (1988). 
Our 0/H abundances at r=0.4R25 versus R///R 25 for 
Virgo spiral galaxies are shown in Fig.^ with circles. A 
trend of decreasing characteristic O/H value with increas¬ 
ing Rif/R 25 is evident for these galaxies. We also derive 
the Hi to optical radius ratio for the field spiral galax¬ 
ies for which 21 cm observations are available. We find 
Rh/I 125 = 1-91 for NGC 628, 1.35 for NGC 1232, 1.14 for 
NGC 2903, 1.02 for NGC 5194, 3.76 for NGC 5236 , 1.75 
for NGC 5457 and 1.49 for NGC 6946, with 21 cm data 
from Shostak & van der Kruit (1984), van Zee & Bryant 
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(1999), Wevers et al. (1986), Shane (1975), Huchtmeier 
& Bohnenstengel (1981), Boulanger & Viallefond (1992). 
Our 0/H abundances at r=0.4R25 versus R///R 25 for field 
spiral galaxies are shown in Fig.^with pluses. Figj^ shows 
that some field spirals have a small Hi to optical radius 
ratios, similar to that in H i deficient Virgo core spirals. 
The field galaxies from our limited sample, with the ex¬ 
ception of NGC 5236, show a similar trend of decreasing 
characteristic 0/H value with increasing R///R 25 as the 
Virgo spiral galaxies. 

Thus, the consideration of the 0/Hp oxygen abun¬ 
dances in HII regions in the Virgo and field spiral galaxies 
suggests that in general the Virgo spirals are indistinguish¬ 
able from the field spirals, both the Virgo periphery and 
core spirals have counterparts among field spirals. There 
are no Hi deficient Virgo spirals with low oxygen abun¬ 
dances. This fact together with the larger number of H i 
deficient galaxies in the core of the cluster causes the ap¬ 
parent segregation of high abundances in the center galax¬ 
ies. 


2.2. Nitrogen 

The N/0 abundance ratio in Hii regions have been deter¬ 
mined using algorithms following below. It was adopted 


N _N+ 

o~o+ 


( 2 ) 


The value is derived from the expression (Pagel et al. 
1992) 


l0g(iV+/0+) = log V//1A6548+A6584 ^ ^ ^ 

VO//]A3726-l-A3729 ^2 

+0.40logt 2 — log(l-I-1.35a;) (3) 


where 

X = 10~‘^net2^^^, (4) 

and He is the electron density in cm“^, t 2 = t[Ar//] is the 
electron temperature in units of lO'^K. It was adopted ^2 
= The t 2 value can be determined from 

the following equations: 

12 + log(0+ + /i7+) = log ^[O^^dA4959+A5007 

Ihp 
1 2'Sl 

—--0.551ogt3, (5) 


12 + \0g{0+/H+) = log ^[0^^1^3726-^A3729 g 

iHp 

0.40logt 2 +log(l +1.35a:), (6) 


O _0+ 0++ 


(7) 

( 8 ) 


where = t[oiii] is in units of lO'^K. Eqs.(||) - (Q) were 
taken from Pagel et al. (1992), Eq.(||) was taken from 
Garnett (1992). Using the value of 0/H derived from 
Eq. (0) and measured line intensities Eqs. (i- ® can be 
solved for The value of rie was adopted to be equal 

to 100 cm“^ for every H ii region. 

The value of t^oii] c an be also found from the following 
expression for tp = ts ( Pilyugin 2001a ) 


i?23 + 3.09 + 7.05P + 2.87P2 
9.90 + 11.86P + 7.05P2 _ 0.583R23 


(9) 


and Eq.(||). The value of ^2 determined from Eqs.( 8 ),(i) 
coincides with the value of ^2 derived from Eqs.(|^),(5)-(B) 
for majority of H ii regions considered. However, for H ii 
regions, in which the bulk of oxygen is in the O'*" stage, 
the two values of t 2 are not in agreement. Eqs. ©Xl)-® 
give a more realistic values of t 2 for these Hii regions. 

Using Eqs. the N/0 abundance ratios were 

derived for all the Hii regions in Virgo and field spi¬ 
ral galaxies for which the [Nil]AA6548,6584 line mea¬ 
surements are available in the works cited in Table 0 - If 
only the [NH]A6584 line measurement was available then 
the total [Nil] AA6548, 6584 line intensity was derived as 


I 


[Af//]AA6548,6584 


= 1.3XI| 


[Ar/7]A6584- 


Fig.|| shows the N/0 versus 0/H diagram for Virgo 
and field spiral galaxies: Virgo core galaxies (NGG 4501, 
NGG 4571, and NGC 4689) are presented as squares, 
in panel a, the Virgo intermediate galaxies (NGC 4254, 
NGC 4321, and NGC 4654) are shown as triangles in 
panel b, and the Virgo periphery galaxies (NGC 4303, 
NGC 4651, and NGC 4713) are presented as circles in 
panel c. Points in all panels represent the field spiral 
galaxies (NGC 628, NGC 1232, NGC 2903, NGC 3184, 
NGC 3351, NGC 5194, NGC 5236, NGC 5457, and 
NGC 6946). 


Panel a of Fig.^ shows that the Virgo core spirals oc¬ 
cupy the lower part of the band outlined by the field spi¬ 
rals in the N/0 - O/H diagram. It looks as if the position 
of Virgo core spirals is slightly shifted (by around 0.05 
+ 0.10 dex) towards higher 0/H values or towards lower 
N/0 values relative to the overall position of field spirals. 
Is this small offset of Virgo core spirals in the N/0 ~ 0/H 
diagram real ? Does it reflect a (marginal) difference be¬ 
tween the evolution of the Virgo core and field spirals? 
It is more likely that this apparent small shift of posi¬ 
tions of the Virgo core spirals in the N/0 - O/H diagram 
is an artifact caused by error propagation in the abun¬ 
dance determination for Hii regions in the core spirals, 
produced by the uncertainties in line measurements (only 
upper limit estimations of [OHI]A5007 line intensities are 
available for those Hii regions), or due to the uncertainty 
in the extrapolation of the P-method into an area of ex¬ 
treme values of the parameters P and R 23 (the value of 
P is equal or close to zero, Fig.0, for all the Hii regions 
of the Virgo core spirals, and the R 23 value is usually less 
than 1 while the calibration is based on the Hii regions 
with R 23 > 2). It should be particularly emphasized that 


t2 = 0.7t3 + 0.3. 
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Fig. 5. N/0 versus 0/H diagram for Virgo and field 

spiral galaxies, a. Virgo center galaxies: NGC 4501, 
NGC 4571, and NGC 4689 (squares) and field spiral 
galaxies: NGC 628, NGC 1232, NGC 2903, NGC 3184, 
NGC 3351, NGC 5194, NGC 5236, NGC 5457, and 
NGC 6946 (points), b. Virgo intermediate galaxies: 
NGC 4254, NGC 4321, and NGC 4654 (triangles) and 
field spiral galaxies (points), c. Virgo periphery galaxies: 
NGC 4303, NGC 4651, and NGC 4713 (circles) and field 
spiral galaxies (points). 


we cannot discard the possibility that this apparent off¬ 
set has an artificial origin, since it is small (around 0.05 
d- 0.10 dex), comparable to the expected errors in the 
oxygen and nitrogen abundances. Therefore error statis¬ 
tics prevent us from any possible data overinterpretation; 
though we stress that more and better data are needed to 
disentangle the actual nature of the offset revealed here. 

Examination of Fig.|| (panels b and c) shows that the 
Virgo periphery and intermediate spiral galaxies occupy 


Table 2. Galaxy Characteristics 


Galaxy 

Name 

Type 

Arm 

Class 

R-efF 

(arcsec) 

Vslrot ,mEix 

(km/s) 

NGC 4303 

4 

9 

62.6 

150 

NGC 4321 

4 

12 

111.5 

270 

NGC 4501 

3 

9 

77.1 

300 


the same band in the N/0 - 0/H diagram as do field 
spiral galaxies, i.e. they appear indistinguishable from field 
spirals in the N/0 - 0/H diagram. 

Thus, the comparison of N/0 abundance ratios in Hii 
regions of the nine Virgo and nine field spiral galaxies 
confirms (or at least is not in conflict with) the above 
conclusion that Virgo spirals result indistinguishable from 
field spirals. If there is an actual difference between abun¬ 
dance properties of the Virgo and field spirals this differ¬ 
ence should be small and appear entirely masked by the 
errors. 


3. Virgo Spirals and Their Multiphase Chemical 

Evolution Models 

Shields et al. (1991) and Skillman et al. (1996) argued that 
the high oxygen abundances in the HI deficient Virgo spi¬ 
rals can be explained by the fact that the infall of metal- 
poor gas onto the galaxy is inhibited in the cluster environ¬ 
ment, and the Virgo core spirals may have evolved more 
nearly in the manner of the closed box ” simple model” of 
chemical evolution. Field galaxies and galaxies on the pe¬ 
riphery of the cluster would have their oxygen abundances 
depressed by infall of metal-poor gas. In the light of our 
results it can be concluded that some field spirals have 
also evolved in the same way as the Virgo core spirals do. 

Another important result of the present study is the 
rather low values of the central oxygen abundances in 
spiral galaxies: the maximum central oxygen abundance 
is ~ 12-|-log(0/H) = 9, about 0.3 dex higher than the 
solar value 12-|-log(0/H)Q = 8.736 (Holweger 2001) or 
12 -|-log(O/H)0 = 8.69 (Prieto, Lambert & Asplund 2001). 
The multiphase multizone models shown by Ferrini et al. 
(1992) predict a maximum value of oxygen abundance 
12 -I- log{0/H) 9.10 dex. This is indicative of that the 
oxygen production used in these models is closer to reality. 
Therefore, one can expect that these models can be used 
to test the influence of the environment, through the effect 
of the infall of gas on the disks of the Virgo galaxies and 
its dependence on distance to the center of the cluster. 

3.1. Multiphase Model Description 


The multiphase mode l used here was fir st applied to 
the Solar Neighborhood ([Ferrini et al. 1992), and to the 
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whole Galactic (MWG) disk (Ferrini et al. 1994 - here¬ 
after FMPD) and bulge (Molla & Ferrini 1995). The same 
model was t hen applied to dis ks ( Molla et al. 1996| , 1999) 
and bulges ( Molla et al. 200^ of some spiral galaxies of 
different morphological types with reasonable success. We 
now present the application of this multiphase chemical 
evolution to three galaxies in Virgo, NGG 4501, close to 
the center, an intermediate galaxy as NGG 4321 and a HI 
normal galaxy NGG 4303, which is in the periphery of the 
cluster, as typical examples for these different zones. 

The data used in our models for these three galaxies 
are presented in Table ^ where we display their morpho¬ 
logical types (T) and Arm Glass in columns (2) and (3) as 
derived by Tully (1988) and Biviano (1991), respectively. 
Effective radii from Tully (1988) are in column (4), with 
the exception of NGG 4303 for which the effective radii 
is from Henry et al. (1992). The adopted distance is 16.8 
Mpc. 

In the multiphase model a protogalaxy is assumed 
to be a spheroid composed of primordial gas with to¬ 
tal mass M(R). For our galaxies the masses are calcu¬ 
lated from the rotation curves (Table |^, column 5) ob¬ 
tained either via radio ( Guhathakurta et al. 1988|) or op¬ 
tical ( Distefano et al. 1990| ) observations. Each galaxy is 
divided into concentric cylindrical regions 1 kpc wide. The 
model calculates the time evolution of the halo and disk 
components belonging to each cylindrical region. The halo 
gas falls into the galactic plane to form the disk, which is a 
secondary structure in the multiphase framework, formed 
by the gravitational accumulation of this gas at a rate 
which depends on the collapse time scale TcoU- 

In our framework an infall dependent on galactocentric 
distance and a two-steps star formation law, by forming 
molecular clouds and then stars, are assumed. 

The infall rate / is inversely proportional to the col¬ 
lapse time scale TcoU , which is assumed to depend on galac¬ 
tocentric radius, through an exponential function with a 
scale length 1: 

TcoiiiR) = ( 10 ) 

The characteristic time scale tq is defined for every 
galaxy as being appropriate to a region equivalent to the 
Solar Neighbourhood in MWG located at a radius Rq 
which is calculated from Reff. This characteristic time 
scale is calculated for every galaxy via the total mass of the 
galaxy. We compute, from the rotation curves, the total 
mass for each galaxy. This defines the collapse time scale 
for each one of them through the expression tq oc Mg T 
( pallagher et al. 1984 ), where Mg is the total mass of the 
galaxy in 10^Mq and T is its age, assumed 13 Gyr. Thus, 
we obtain tq for each galaxy, through the proportion: 
To = TQ{Mg,gal/Mg^MWG)~^^'^, the value of r© = 4 Gyr 
having been determined in Ferrini et al. (1992) for the 
Solar Region. 

Following the relation between collapse time scale and 
mass, it is evident that TcoU must be variable with the 
galactocentric radius. If the total mass surface density fol¬ 
lows an exponential law, such as the surface brightness 


Table 3. Gomputed Models 


Galaxy 

Name 

Model 

Ro 

(kpc) 

To 

(Gyr) 


en 

NGG 4303 

B 

6 

8 

0.22 

0.008 

NGG 4303 

M 

6 

16 

0.20 

0.005 

NGG 4321 

B 

10 

7.5 

0.45 

0.160 

NGG 4321 

M 

10 

15 

0.45 

0.320 

NGG 4501 

B 

8 

6.5 

0.40 

0.150 


does for spiral disks, the resulting TcoU must increase along 
the radius with a scale length I ex Re (where Re is the 
scale-length for the surface brightness radial distribution). 
We use I = 4 kpc but we must take into account that Re 
usually decreases for later type galaxies and is larger for 
the earlier ones, and therefore the selection of this param¬ 
eter have some uncertainties. 

In the various regions of the galaxy we allow for differ¬ 
ent phases of matter aggregation: diffuse gas (g), clouds 
(c, only in the disk), low-mass (si, m < dM©) and massive 
stars {s 2 ,m > 4M©), and remnants. The mass in the dif¬ 
ferent phases changes through several conversion process: 

1. Star formation by the gas spontaneous fragmentation 
in the halo (oc Kgjj^) 

2. Formation of the disk by the accumulation of diffuse 
gas from the halo (oc fgn) 

3. Gloud formation in the disk from diffuse gas (oc gg^'^) 

4. Star formation from cloud-cloud collisions (oc Hc^) 

5. Induced star formation via massive star-cloud interac¬ 
tions (oc acS' 2 ) 

6 . Diffuse gas restitution from these two star formation 
processes 

The rates for these processes are proportional to pa¬ 
rameters /i, H, a and K, which depend on galactocentric 
radius through the equations computed in FMPD. The 
proportionality factors of those equations are the corre¬ 
sponding efficiencies of processes, that is the probability 
of cloud formation, e^, of cloud—cloud collision, Cij, and of 
the interaction of massive stars with clouds, ta in the disk, 
and the efficiency to form stars in the halo, ck- For each 
galaxy these characteristic efficiencies must be chosen, al¬ 
though the term of the induced star formation is associ¬ 
ated to a local process and, as a result, its coefficient is 
considered independent of both position and morphologi¬ 
cal type. The last term ex is also assumed constant for all 
halos, thus being independent of morphological type. The 
other rates e^ and en, depend however on the Hubble type 
and/or the arm class, and their adopted variation range 
is determined following the arguments given by Ferrini & 
Galli (1988), as discussed in Molla et al. (1996, 1999). 


In order to reduce to a minimum the freedom degree 
due to the high number of parameters, the model has been 
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applied following a precise strategy: first, it was used for 
the Solar Neighbourhood, followed by the MWG model, 
to determine the characteristic parameters of MWG. Then 
we afforded the final generalization for some other spiral 
galaxy, by assuming that the two efficiencies Cq and ck do 
not change from galaxy to galaxy. Therefore, we allow to 
change from galaxy to galaxy, only the characteristics col¬ 
lapse time scale, tq, depending on the total mass and the 
two efficiencies, and ch, following the morphological 
type. The selection of each one of these input parame¬ 
ters in our models will be the matter of discussion in the 
following subsection. 

The adopted initial mass function (IMF) is taken from 
Ferrini et al. (1990). The enriched material proceeds from 
the restitution from dying stars, considering their nucle¬ 
osynthesis, their IMF (and the delayed restitution) and 
their final fate, via a quiet evolution, or Type I or II su¬ 
pernova (SN) explosions. Nucleosynthesis yields are from 
Renzini & Voli (1981), and Woosley & Weaver (1995) 
for low-mass and intermediate, and massive stars, respec¬ 
tively. The type I supernova explosions release mostly iron 
following Nomoto et al. (1984) and Branch & Nomoto 
(1986) at a slower rate, by implying that the iron appear 
at least 1 Gyr later than the so-called a-elements (Oxygen, 
Magnesium ...) ejected by the massive stars. 

3.2. Computed models 

The collapse time scales have been chosen according 
to the above paragraph and their values are shown in 
Table For the efficiencies, they usually turn out to be 
larger for earlier morphological types and lower for the 
later ones in MFD96, and this is useful for a first se¬ 
lection of the values for a given galaxy. The fine tun¬ 
ing is then performed, once the collapse time scale pa¬ 
rameters To and I chosen, when and en are con¬ 
strained to reproduce the observational data. We use 
for this task the observed radial distributions for atomic 


gas Hi 

( 

Warmels 198S, Cayatte et al. 1990]), molecular 

gas H 2 , 

( 

Kenney & Young 1988], 

Kenney & Young 1989, 

Rownd & Young 1999), star formation rate surface den- 


sities, obtained from Ha fluxes, from Kenney & Young 
(1989) and Rownd & Young (1999), and oxygen abun¬ 
dances of the gas phase 12 -I- log{0/H) derived above. 
For NGG 4321 we also used the most recent data on Ha 
fluxes and H i and H 2 radial distributions from Knapen & 
Beckman (1996) and Knapen et al. (1996). 

The set of models is shown in Table ^ where we give, 
for the region located at a distance Rq, column ( 6 ), equiv¬ 
alent to the Solar Region, the corresponding values of the 
characteristic collapse time scale, column (7) and the two 
efficiencies and en, columns ( 8 ) and (9). 

3.3. The present-day radial distributions 

In the following graphs we show the results derived from 
the chemical evolution models defined by the input data 



t (Gyr) 


Fig. 9. Infall rate evolution for two regions in each galaxy, 
inner (4 kpc) -strong line- and outer (14 kpc) -slight line- 


of Table |3|, together with the corresponding observational 
data. FigTg, 0 and ^ show for each of the studied galaxies, 
the radial distribution of the oxygen abundance in panel 
(a); the star formation rate surface density normalized to 
its value for each equivalent solar region 4>o in (b), the 
gas surface density of diffuse Hi gas in (c), and the gas 
surface density of molecular gas H 2 in (d). 

Abundances are really well reproduced for the three 
galaxies. By using the results derived in this work, the fit 
improves extraordinarily, in comparison with models from 
Molla et al. (1999), where abundance estimates from the 
R 23 -method were used, thus supporting our new oxygen 
abundances. 

Atomic gas radial distributions are quite well repro¬ 
duced for all three galaxies, although they have a high gas 
surface density in the outer disks in comparison with the 
data. This difference between data and models does not 
appear in NGC 4303 models. For Model B of NGG 4321 it 
is ~ 1 — 1.5Mq/pc^ for R > lOkpc, while for Model M it 
does not exist. This implies that this last model is better 
able reproduce the outer disk, while the first one better fits 
the data in the inner regions. For the NGG 4501 model, the 
difference between data and models appears at i? > 8kpc 
and takes values comparable to the mean density in this 
outer disk. Summarizing, the discrepancy between models 
and data increases with decreasing distances to the center 
of the cluster, simulating perfectly the behavoir of the de¬ 
ficiency of the gas. This observed depletion of the diffuse 
gas from the outer regions of disks needs some mechanism 
producing this gas loss, which is not included in our mod¬ 
els. Therefore we cannot reproduce this special feature of 
Hi deficiency in the Virgo outer disks. 
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Fig. 10. Relation of oxygen abundances with the gas frac¬ 
tion fj,. The solid line represent the closed box model pre¬ 
diction, while points correspond to the radial distributions 
of oxygen abundances for the three galaxies at the present 
time as indicated in the figure 


Molecular gas radial distributions differ slightly in the 
inner regions, but always close to the error bars. The 
uncertainties in the estimation of the molecular cloud 
masses, which are obtained from CO intensities through 
a calibration parameter y, are large because the de¬ 
pendence of X on. metallicity or galactocentric distance 
( Verter fc Hodge 1995 , Wilson 1995 ) in each galaxy is not 
completely well determinated. 

These good results in the predicted gas radial distri¬ 
butions are, in fact, the consequence of searching mod¬ 
els fitting these two distributions. Therefore the largest 
differences for a galaxy may only appear in the oxygen 
abundance and star formation rate radial distributions. 

As we have already pointed out, abundances show a 
good behaviour, consistent with the data, while the shape 
of the star formation rate surface density radial distribu¬ 
tions are well fit, showing exponential shapes, with the 
exception of Model B for NGC 4321 (panel d. Fig. 
This model B predicts a star formation rate lower than 
the data in the inner regions. This also appears in Fig. || 
where a maximum rate is produced at 3-4 kpc and it im¬ 
plies that the model has a strong evolution, maybe more 
rapid than the actual evolution. 

Thus, it seems that, following Figs. |[ our models 
reproduce the abundance data in Virgo galaxies well, and 
therefore we can use these models to explain the time evo¬ 
lution of Virgo galaxies. In order to see clearly the effects 
of the time evolution we show in Fig.| the infall rate for 
two regions, located at 4 and 14 kpc, and for model B of 
the three galaxies. We see that the infall rate is clearly 


higher in the inner region of NGC 4501 compared to the 
two other galaxies, NGC 4321 rests in middle of the two 
others, while NGC 4303 has a slow infall even in this cen¬ 
tral region. For the outer regions, the infall rate continues 
to be larger for NGC 4501, intermediate for NGC 4321 
and low for NGC 4303, following the same trend as in the 

4 kpc radial zone, but differences among galaxies are now 
smaller. All zones and galaxies show a similar and low 
infall rate at the present time. The most important dif¬ 
ference resides in the fact that the infall rate had a very 
strong maximum during the first two Gyr for the inner 
diks in both NGC 4501 and NGC 4321, which finished 
at ~ 5 Gyr, while NGC 4303 maintains a slow evolution, 
with an infall rate higher than that of the two others from 

5 Gyr up to now. 

These different infall rate histories have important ef¬ 
fects on the oxygen abundances obtained along the galac¬ 
tic disk for these galaxies. We show in Fig.|l^ the oxygen 
abundances, as 12 -|- log (0/H), for each radial region and 
galaxy vs the gas fraction /r = Mgas/Mdisk- In this graph 
the solid line represents the closed box model prediction, 
while points correspond to the three disk radial distri¬ 
butions. The positions of NGC4303 (circles), NGC4321 
(triangles), and NGC4501 (squares) in the ^ - 0/H dia¬ 
gram show that these galaxies evolve with some gas infall. 
We may see that the outer regions of NGC 4303 are close 
to the closed box model line, simulating a galaxy with¬ 
out mass exchange with the environment, while the inner 
regions have abundances smaller than the ones expected 
from their gas fraction according to the closed box model. 
NGC 4321 and NGC 4501 show a strong decrease in its 
abundances for their inner regions, and there is, even for 
the outer ones, some dilution effect. All these regions have 
very low gas fractions, which according to the closed box 
model should correspond to the highest abundances, 9.0- 
9.10; however this is not the case, and the abundances are 
lower than these values. 

Thus, our models show that both Hi normal and Hi 
deficient galaxies have their oxygen abundances depressed 
by infall of metal-poor gas. Indded, the model for the 
NGC 4501 (a representative of the Hi deficient Virgo core 
spirals) shows that the inner regions of NGC 4501 have 
very low gas fractions and therefore even a slow infall rate 
at the present epoch results in a dilution effect, and the 
evolution of Virgo core spirals are somehow similar to that 
of field spirals. Taking into account that our models repro¬ 
duce the data well for the three galaxies, in particular the 
gas densities and abundances, this conclusion does not 
depend on the goodness of the hypotheses or theoretical 
assumptions. 

4. Discussion 

In order to make sure that our models are in agreement 
with the whole set of new results for the cluster galaxies 
and in particular for Virgo, we will summarize these ob¬ 
servations. Solanes et al. (2001) have derived the atomic 
content for a total of 1900 galaxies in 18 nearby clus- 
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ters, aiming to study possible connections between the 
gas deficiency and the characteristics of the galaxies or 
their environment. With these data it is well established 
that the proportion of gas-poor galaxies increases contin¬ 
uously toward the center of the cluster for a large num¬ 
ber of them. Furthermore, there are clear indications that 
the removal of H i occurs mostly in the outer regions of 
disks, thus implying that the Hi disk sizes are reduced 
(Warmels 1988, Cayatte et al. 199C, Cayatte et al. 1994), 
and that the degree of the H i depletion is related to the 
morphological type, early types and dwarfs being more 
easily emptied of gas than intermediate types. 

On the other hand, the analysis of the spectroscopic 
catalog of distant clusters performed by Dressier et al. 
(1999) and Poggianti et al. (1999) show that the star for¬ 
mation rate has been quenched, rather than enhanced, 
in these galaxies. This conclusion has been also reached 
by Caldwell et al. (1999), who found that star forma¬ 
tion takes place in the inner regions of these galaxies, 
but it is suppressed in the outer ones. Moreover the col- 
ors of these galaxies ([Gavazzi, Boselli, &: Kennicutt 19^ 


Gavazzi et al. 1998) do not differ significantly from those 
of field galaxies for the same morphological types, and the 
molecular gas content, usually more concentrated in the 
inner disks, show normal masses (Kenney & Young 1988, 
1989). 

Thus, the scenario necessary to reproduce all observa¬ 
tions must be able to eliminate the gas of the outer regions, 
but not in the inner disk, and at the same time, it has to 
form stars at the early times of evolution, in order to ob¬ 
tain the red colors observed and the star formation rates 
estimated for the present time. The most extended idea 
is that the mechanism producing this removal of gas is 
the ram pressure stripping due to the orbit of each galaxy 
around the core of the cluster. Under this assumption the 
low density diffuse gas may be eliminated when the galaxy 
passes close to this core, while the higher density molecu¬ 
lar clouds located in the inner disk should not be affected. 

Recently, Vollmer et al. (2001) have performed some 
N-body simulations in order to study the effects of this 
mechanism on a galaxy suffering ram pressure stripping 
by its motion within the cluster. They find that the move¬ 
ment of the galaxy around the cluster center provokes an 
acceleration of the clouds located on the surface of the 
disk, which, in turn, produces a high surface density in 
the center of the galaxy. Thus, the star formation can suf¬ 
fer an enhancement in the inner disk. After this phase, the 
gas of the outer regions is removed and, without this fresh 
gas reservoir, the star formation decreases. 

Within this scenario, the high density produced in the 
inner regions of the galaxy is the consequence of a strong 
infall rate on a short time scale. This is exactly what hap¬ 
pens in the NGC 4501 model, where the inner disk suf¬ 
fers a large infall rate in a short time, by creating stars 
very rapidly. The consequence is the consumption of the 
gas, and therefore star formation stops. For a Virgo pe¬ 
riphery galaxy, as NGC 4303, the infall is much slower, 
maintaining a more constant star formation rate along 


its evolution, similar to the typical one of a field galaxy. 
The model for NGC 4321 has a behavior intermediate be¬ 
tween both of them. This explains why the star forma¬ 
tion rate is, for the present time, higher for NGC 4303, 
which still has enough gas, while this rate has decreased 
for NGC 4321 from its maximum value, and it is very de¬ 
pressed for NGC 4501, which has almost no gas. These 
model predictions are supported by the observations from 
Koopmann & Kenney (1998) and Rownd & Young (1999) 
who find that the star formation efficiency decreases with 
increasing Hi deficiency among the Virgo cluster galaxies. 
In fact the Ha luminosities give values logL = 8.38, 7.97 
and 7.96 Lq for NGC 4303, 4321 and 4501, respectively, 
following these latter authors. 

We would like to add some other considerations about 
our models. Koopmann & Kenney (1998) found that the 
morphological classification for the Virgo cluster galax¬ 
ies does not follow the same scheme as the one used for 
the field galaxies. In this last case the concentration in¬ 
dex (measured as the bulge to disk luminosity ratio) is 
correlated with the star formation in the arms, while in 
the cluster galaxies it does not occur: due to the Hi defi¬ 
ciency in the outer disk, the concentration index increases 
for a given Ha flux, as a measure of the star formation 
rate. Thus if it corresponds to a Sc galaxy, the concentra¬ 
tion index would define the galaxy as Sa. This fact implies 
that the Virgo galaxies would be classified as earlier types 
whereas their star formation rates would correspond to 
later types of galaxies. In this sense, it appears remarkable 
that our models reproduced the observed characteristics 
of the three selected galaxies with efficiency values lower 
than those typically used for their morphological types, 
and closer to the values corresponding to one unity later 
morphological types. 

5. Conclusions 

The oxygen and nitrogen abundances in H ii regions in the 
nine Virgo spirals of the sample from 
and in nine selected field spiral 
determined with the recently suggested P - method. 

It has been confirmed that there is an abundance seg¬ 
regation in the sample of Virgo spirals in the sense that 
the HI deficient Virgo spirals near the core of the cluster 
have higher oxygen abundances in comparison to the spi¬ 
rals at the periphery of the Virgo cluster. In general, both 
Virgo periphery and Virgo core spirals have counterparts 
among field spirals. Some field spirals have a small Hi to 
optical radius ratios, similar to the one in Hi deficient 
Virgo core spirals. It has been concluded that if there is 
an actual difference in abundance properties between the 
Virgo and field spirals this difference should be small or 
entirely masked by the errors. 

These abundance results have been analyzed with the 
multiphase chemical evolution model, applied to three 
galaxies NGC 4501, 4321 and 4303, considered as typical 
examples of galaxies located at the center, at intermedi¬ 
ate locations, and at the periphery of the cluster. These 


Skillman et al. 1996 
galaxies were re- 
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models are shown to be able to reproduce the observa¬ 
tions of abundances, star formation rate and diffuse and 
molecular gas densities. Following these models the infall 
rate was strong at early times and it is very low now for 
the Virgo core galaxies, but with features continuously 
smoother for the periphery galaxies. The consequence of 
the infall of the gas is the dilution of the elements, and 
therefore lower abundances than those expected with the 
closed box model, for the gas fractions in the center galax¬ 
ies, are predicted. Thus, the final abundances turn out to 
be similar to the ones present in field galaxies. 

The scenario of the ram pressure stripping usually 
claimed to produce the observed characteristics in the 
cluster galaxies, such as a H i deficiency, which is larger for 
the galaxies closer to t he cluster, was rece ntly simulated 
with N-body models by Vollmer et al. 2001 . Following this 
work, a large infall rate is produced in the inner regions of 
the disk of these galaxies, as a consequence of their move¬ 
ment towards the cluster core, with a posterior phase of 
quenching of the star formation; these results appear con¬ 
sistent with our conclusions. 
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Fig. 1. Oxygen abundances and values of excitation parameter P versus galactocentric distance for Virgo spiral 
galaxies. The oxygen abundances determined with the R 23 - method calibrated by Zaritsky et al. 1994 are shown by 
pluses, the best fits t o these data ar e presented by dashed lines. The oxygen abundances determined with the P - 
method calibrated by Pilyugin 2001a are shown by circles, the best fits to these data are presented by solid lines. The 
values of excitation parameter P are shown by triangles. The galactocentric distances are normalized to the isophotal 
radius. 
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Fig. 2. The same than Fig.^ for the field spiral galaxies. 
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Fig. 6. Predicted Radial distributions at the present epoch for NGC 4303: (a) Oxygen abundance 12 + log{0/H)-, (b) 
Star Formation Rate surface density normalized to the equivalent to solar region value4'/'I'o; (c) diffuse or atomic gas 
Hi and (d) molecular gas H 2 surface density. Model is represented by the lines while symbols are the observational 
data 
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Fig. 7. Predicted Radial distributions at present time for NGC 4321 with the same panels and symbols than in Fig. ^ 
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Fig. 8 . Predicted Radial distributions at present time for NGC 4501 with the same panels and symbols than in Fig. ^ 













